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We present an investigation of the effects of constant but random shifts of the ground hyperfine 
qubit states in the setting of quantum computing with ion doped crystals. Complex hyperbolic 
secant pulses can be used to transfer ions reliably to electronically excited states, and a perturbative 
approach is used to analyse the effect of ground state hyperfine shifts. This analysis shows that 
the errors due to the hyperfine shift are dynamically supressed during gate operation, a fact we 
attribute to the AC Stark shift. Furthermore we present an implementation of a controlled phase 
gate which is resilient to the effects of the hyperfine shift. Decoherence and decay effects are included 
in simulations in order to show that a demonstration of quantum gates is feasible over the relevant 
range of system parameters. 

PACS numbers: 03.67.Lx, 32.80.Qk, 03.65.Yz 



I. INTRODUCTION 

An interesting class of proposals for quantum comput- 
ing is provided by solid state devices with dopants as the 
active medium [H, 0, 0, 0- In this article we focus 
on the Rare Earth Quantum Computing (REQC) sys- 
tem proposed by Ohlsson et al. 4]. In this proposal 
qubit states are encoded in hyperfine levels of rare earth 
ions doped in an inorganic crystal host. The qubit lev- 
els are coupled via an excited state accessible by optical 
excitation. Due to the large inhomogeneous broadening 
of the excited state, qubits may be selected in frequency 
space and identified by the magnitude of the inhomoge- 
neous shift. In practice qubits are prepared as narrow 
structures in a hole burnt structure within the inhomo- 
geneous profile and readout is achieved using absorption 
spectroscopy (see Fig. [IJ. Each qubit consists of an en- 
semble of ions centered at a certain transition frequency. 
For this reason operation of the REQC system calls for 
pulses that are robust to variations in detuning inherent 
in the proposal. Qubit interaction is mediated by the 
dipole-dipole coupling. The rare earth ions within the 
crystal acquire a change in static dipole moment when 
in the excited state |e). This causes the excited state of 
the surrounding ions to shift in energy, and provided the 
shift exceeds a certain threshold, resonant driving of the 
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FIG. 1: Left: simplified level diagram of a single ion. Right: 
qubits are prepared as antiholes in a hole burnt structure. 
To ensure appropriate interactions between qubits the system 
must first be properly initialised Q. 
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FIG. 2: Implementation of a controlled NOT gate. Pulses 2-4 
implement the NOT gate on the target qubit. If the control 
qubit is initially in the |0) state it is promoted by pulse 1 to 
the excited state and pulses 2-4 become off resonant. If, on 
the other hand, the control starts out in the |1) state pulse 1 
is off resonant and the NOT gate on the target ion proceeds. 



\i) <-* |e) transition is blocked. This effect can be used 
to implement a fully entangling gate |4| as illustrated in 
Fig.H 

Previous theoretical works have dealt with the problem 
of transferring population reliably in spite of the finite 
width of the frequency channels [5| and the use of phase 
compensating pulses to eliminate the effects of unknown 
detuning dependent phases acquired during pulses |6( . In 
the laboratory qubit initialisation ,7] and basic manipu- 
lations have been demonstrated with good fidelity and 
it appears that a fully entangling gate will soon be within 
reach. 

In this paper we focus on an effect that has previously 
not been addressed in the REQC system. Apart from the 
excited state shift due to the finite qubit channel width 
which is chosen in the initialisation procedure, there is 
also a broadening of the hyperfine qubit levels with a 
width of about 30 kHz. The constant but random shift of 
the hyperfine levels within this range leads to dephasing 
between the |0) and |1) states and is totally destructive to 
previously proposed gate implementations. In this paper 
we propose gate sequences that correct for the ill effects 
of the hyperfine shift. 
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The detrimental effect of atomic decay is often ne- 
glected in theoretical investigations of REQC. Although 
this may be well justified for short pulse sequences with 
negligible integrated population of the excited state, it 
becomes important in longer error compensated schemes. 
We thus present simulations which include a description 
of the main decoherence mechanisms. 

The paper is organised as follows. In Section [H] we 
give a theoretical description of the complex hyperbolic 
secant pulse and discuss the effects of the hyperfine shift. 
In Section IIIII we investigate the effects of the hyperfine 
shift on a previously reported method for performing ar- 
bitrary qubit rotations on single qubits. In Section lTvl we 
propose a phase compensated implementation of a con- 
trolled phase gate which is robust to the effects of the 
finite channel width as well as the hyperfine shift. Sec- 
tion deals with the main contributions to decoherence 
and we perform a full simulation of our gate sequence 
taking these effects into account. 



II. THE COMPLEX HYPERBOLIC SECANT 
PULSE 

The complex hyperbolic secant pulse offers population 
inversion over a broad range of dctunings 9] and is there- 
fore perfectly suited for rare earth quantum computing. 
A detailed theoretical analysis given in ^J] uses the Bloch 
sphere picture of two-level dynamics and shows that the 
secant pulse, in fact, satisfies very well the adiabaticity 
criterion during an effective chirp through atomic reso- 
nance. The Bloch sphere approach does not reveal the 
acquired global phase on the two coupled quantum states, 
which is important in the present quantum information 
setting since in REQC we always address a subspace 
TL = {\e),\i)} of the full three level system and hence 
a global phase acquired during the |0) <-> |e) transitions 
results in a relative phase between the |0) and |1) qubit 
states. In the following we present a perturbative treat- 
ment which reveals the global phases acquired during a 
Sech pulse. 

The complex Rabi frequency of the Sech pulse is 

n scch (t) = ^(sech^t)) 1 "^ = Cl R (t)eW*\ (1) 



with the real function 



Q R (t) = Q sech{/3t), 



(2) 



and the instantaneous frequency given by the time 
derivative of the phase 



tp(t) = ^tanh(/3t). 



(3) 



During the pulse the evolution is determined by the time 
dependent Hamiltonian. For a given ion with a transi- 
tion frequency shifted by A with respect to the centre 
frequency of the channel we arrive at the Hamiltonian 



H{t) = (A + 0(t))|e)(e| + \£l R {t){\i)(e\ + he. 



(4) 



with i — 0, 1 and K = 1. Note that we are working in 
an accelerated frame rotating at ip(t). Diagonalising the 
instantaneous Hamiltonian we find the time dependent 
energies 



E± {t) = A + y(t)± v /(A + y(t))2 + ^) ^ 



and the corresponding eigenstates 

2E±(t)\e)+Sl R (t)\i) 



l±> 



(6) 



Expanding the Schrodinger equation in the time depen- 
dent eigenbasis we obtain the equations of motion 



ic± = E ± c± - (±|=i=>c T . 



(7) 



In the adiabatic limit the coupling term (±|=p) tends to 
zero and the action of the Sech pulse becomes an adia- 
batic transfer. The equations of motion are then simply 
solved by integrating the time dependent energies 



4 0) (*) 



-if*E±dt' 



(8) 



We now include the effect of the diabatic correction as 
a perturbation. Switching to the interaction picture we 
have to first order 



Ui = l-i / Vidt' 



(9) 



The interaction picture potential in the basis {| — ), |+}} 



is 



Vi = 



i£ e i Jo E+-E-dt' 



-i£ e ij£ E --E+dt^ 



with 



fl R (A + ip) — ifCl R 
2((A + fY + 0|) 



(10) 



(11) 



In Fig. we plot the Bloch sphere trajectories dur- 
ing a Sech pulse computed to zeroth and first order and 
compare with a numerical integration of the equations 
of motion. We observe good agreement between the nu- 
merical result and the first order perturbative result so 
we can dispense with higher order corrections. Note that 
as long as we are restricting ourselves to the two dimen- 
sional subspace H = {|e), \i)} the effects of the hyperfine 
shift can always be absorbed into the effect of the excited 
state shift. Only when we regard the full three level sys- 
tem do the effects of the hyperfine shift materialise. 



III. ARBITRARY QUBIT ROTATIONS 

The Sech pulse described in the previous section is only 
able to implement effective 7r-rotations. A method to 
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TABLE I: Phases tpo and (pi of the bichromatic Sech-pulse 
fields. The first two pulses apply a phase of 6 to the |1) 
state. Pulses 3-4 compensate the detuning dependent phase 
acquired during the first two pulses. 
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FIG. 3: (Color online) Bloch sphere trajectories during a Sech 
pulse with pulse parameters Qo — 4 MHz, fj, = 3, j3 = 1.28 
MHz and T — 1.5 [is. The North and South poles represent 
the states |e) and \i) respectively. The solid curve is the result 
of a numerical integration with a standard ODE solver. The 
dot-dashed curve (left) represents the adiabatic +) eigenstate 
of Eq. @ and the dashed curve represents the perturbative 
result of Eqs. H- I1H . The inset shows the ODE result in a 
frame rotating at constant frequency. 



perform arbitrary single qubit rotations has been sug- 
gested in By simultaneously applying two fields one 
can selectively address one of the two 'bar' states defined 
by 
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We selectively address the |0) <-> |e) transition by simul- 
taneously applying two fields with the |0)-|e) and |l)-|e) 
optical transition frequencies and with complex Rabi fre- 
quencies firj(i) = ^sech(i) and f2i(i) = e 4 T2s GC h(i) with 
ip = a. In order to address the |1) «-> |e) transition we 
choose the relative phase such that ip = a + %. 

Using Sech pulses we can selectively apply a phase to 
the |1) state, i.e. effect the evolution 



U 



1 

e lt 



(13) 



in the basis {[0), |T)}. The phase e 10 is applied to the 
|T) state by first applying a two-color Sech-pulse to the 
1) «-* |e) transition. A second Sech-pulse is then applied 
with the phase shifted by e relative to the first pulse 
returning the ions to the |T) state with a geometric phase 
e l9 . During the two pulses the |T) state in addition to 
the geometric phase acquires a phase </?(A) dependent 
on the detuning from the qubit channel center. Now, 



FIG. 4: In the logical basis inhomogeneous broadening of the 
qubit levels leads to dephasing(left). In the 'bar' basis the 
effect is a coupling of the |0) and |I) states(right). 



applying two Sech-pulses with a relative opposite phase 
to the |0) <-> |e) transition the |0) state picks up the same 
detuning dependent phase y(A), but no geometric phase, 
and <p(A) may be taken as a global phase and thus be 
disregarded. The phases of the two fields Slo(t) and Oi (t) 
during this four-pulse sequence are summarized in Table 

Q 

In the logical basis the evolution described by Eq. 113|) 
is equivalent to 



U = e ie ' 2 



cos (9/2 ie-w sin 0/2 
ie lv sm6/2 cos6>/2 



(14) 



which is an arbitrary rotation in the qubit space. 



A. Effects of inhomogeneous broadening 

The Sech pulse offers robustness against variations A 
in the resonant optical excitation frequencies. We now 
include in our formalism a term <5|1)(1| accounting for 
the inhomogeneous shift of the hyperfine ground levels. 
The Hamiltonian then becomes 



H(t)= (A + 0(i))|e)(e|-£|l)(l| 
+ifi o (t)|e>(0|+hc. 
+ |fi 1 (t)|e)(l|+hc. 



(15) 



Like A, the shift 8 varies for the different ions, and the 
goal is to implement gates which are also resilient to 
this variation. The Hamiltonian in the basis {|e), |0), |1)} 
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takes the form 

H(t)= (A + 0(t))|e)(e|-|(|O)(O| + |l)(l|) 
+^(n (t)+e la n 1 (t))\e)(Q\ + hc. 

+ 5 ^(noW-e fa n 1 (t))| e >(i| + hc. 

+|(|0)(1| + |1)(0|). (16) 

Since the inhomogeneous shift of the hyperfine ground 
states is of the order 5 ^ 30 kHz while the Rabi fre- 
quencies may be as large as 4 MHz we may treat the 
effect of hyperfine broadening perturbatively. We split 
the Hamiltonian of Eq. l|16f) into two parts 

H(t)=H a (t) + V, (17) 

with 

H = (A + <p(t))\e){e\ 

+ ^(n Q (t)+e la n 1 (t))\e)(d\+hc. 

+ ^(tto(i)-e M fi 1 (t))|e)(l|+hc. (18) 

and 

y=f(-|0)(0|-|I)(l| + |0)(l| + |l)(0|). (19) 

Since we only address one transition at a time the dynam- 
ics due to H is restricted to the subspace H — {|e), \i)} 
and is thus described by the methods developed in Sec- 
tion [H] whereas V causes leakage of population between 
|0) and |1). Resorting to Eqs. © and (fT5 |) — (fTU |) we find 
that when driving for example the |0) <-> |e) transition 
this leakage is described by 

Ur = ~f 1 "* e-*!t'^"dt\ (20) 

with SIr and E + given by Eqs. 10) and (JSJ) respectively. 

The population of the |1) state during the pulse is 
shown in Fig. [SJ Note that with the given parame- 
ters a naive estimate of the final population would be 
P ~ ST = 30 kHz 1.5 /xs = 0.045 which is an order of 
magnitude greater than the observed value P — 0.0043. 
Intuitively we may understand this as an effect of the AC 
Stark shift. Once the field builds up sufficient strength 
the |0) level is shifted in energy and the coupling S is 
too weak to drive the transition between non-degenerate 
states (see Fig. [HJ ■ We thus understand the form of the 
curve in Fig.[5]as follows: at the beginning of the pulse we 
observe ordinary Rabi oscillations between the |0) and |T) 
states. At t = 0.59 [is we find E + = 106 and we enter the 
regime E + ^> 5 where the curve begins to depart from 
the ideal sinusoidal shape. At t = 0.82 /is E+ = 1006 
and the curve flattens out and begins to oscillate due to 
the complex exponential of Eq. I|20|) . 
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FIG. 5: (Color online) Population leakage due to S calcu- 
lated perturbatively(dashed line) and by numerical integra- 
tion(solid line). Parameters are S7o = 4 MHz, A = 100 kHz, 
S = 30 kHz, = 1.28 MHz and (i = 3. 




FIG. 6: When there is no field the weak coupling 8 cannot be 
ignored, but in the presence of a strong field the coupling of 
the hyperfine levels is suppressed by the AC Stark shift. 



IV. A CONTROLLED PHASE GATE 

In the previous section we saw how the AC Stark shift 
suppresses the dynamics due to 5 during single qubit ro- 
tations. Experimental limitations place an upper limit 
on the laser power so the suppression cannot be made 
complete. When working with two-qubit gates we may 
restrict ourselves to a single fully entangling gate and 
thus have a little more freedom in choosing the imple- 
mentation. In the following we shall consider a robust 
implementation of a controlled phase gate. 

A naive implementation of a controlled phase gate us- 
ing the dipole blockade effect is given in Table ITU The 
table is to be read in the following way. First the control 
qubit is rotated by a Sech pulse through n radians about 
an equatorial axis of the {|e) c , |0) c } Bloch sphere. Sec- 



TABLE II: Naive implementation of a controlled phase gate 
with all pulses implemented as single color Sech pulses. Pulses 
2-3 implement an effective 27r-rotation on the target qubit. 
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TABLE III: Phases <po and (pi of the single color fields during 
an implementation of a single qubit NOT gate. 
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ondly the target qubit is rotated by two consecutive Sech 
pulses through 2tt radians about an equatorial axis of 
the {\e) t , |l)t} Bloch sphere and finally the control qubit 
is rotated through tt radians retracing its path on the 
{|e) C! |0) c } Bloch sphere. Given an arbitrary initial state 
\ip) = Coo|00) + Coi|01} + Cio|10) + Cii|ll) the first pulse 
excites the |0z) components thereby blocking the subse- 
quent excitation of the |01) component. The second pulse 
thus only addresses the |11) component, applying a phase 
of —I. The final pulse returns the \0i) components to 
their initial state. 

Two main sources of errors affect the pulse sequence of 
Table HTl First, the finite channel width causes a detun- 
ing dependent phase on each transition. This defect can 
be remedied by applying phase compensating pulses 0. 
Secondly, the hyperfine broadening causes dephasing of 
the qubits during gate operation. It is well known from 
nuclear magnetic resonance studies that such errors can 
be suppressed by rapidly flipping the two levels. This 
approach is not suitable for REQC since we do not have 
direct access to the |0) <-+ |1) transition. Furthermore fast 
operations are unsuitable since they do not fulfill the adi- 
abaticity criterion. We can, however, perform a robust 
bit flip (or NOT gate) using the three-pulse sequence 
shown in Table ITTT1 4]. Since the Sech pulse must fulfill 
the adiabaticity criterion the bit flip is relatively slow. 
Rather than using the standard 'Bang Bang' technique 
we must use the NOT sequence sparingly while retaining 
some degree of refocussing. To this end we propose the 
pulse sequence shown in Table HVI The sequence consists 
of a phase compensated controlled phase gate implemen- 
tation interrupted by a single refocussing pulse sequence. 
The effect of the different pulses is summarised in Fig. 
In order to evaluate the gate performance we use the 
fidelity defined by 



T= \{ll>in\UlU\ll> in )t 



(21) 



where \tpi n ) is the initial state, U the evolution operator 
implemented by the pulse sequence and Uq the ideal op- 
erator implementing the gate. The fidelity may also be 
defined as the average of T or the minimum of T over the 
Hilbert space. We are only interested in a single charac- 
teristic of the gate performance, and we shall evaluate 
Eq. 12I|) for a single superposition input state. 

The fidelity of the full gate implementation as a func- 
tion of the hyperfine shifts of the two qubits is plotted in 
Fig. |SJ We observe good fidelities over the relevant pa- 
rameter range, in particular the fidelity is virtually inde- 
pendent of the shift of the control qubit. In fact we expect 



TABLE IV: Implementation of a robust controlled phase gate. 
Pulses 4 and 5 compensate detuning dependent phase errors 
caused by pulses 2 and 3. Pulses 7-12 as well as 19-24 imple- 
ment a NOT gate on both qubits. Pulses 13-18 compensate 
detuning dependent errors caused by pulses 1-6. Note that 
all pulses are implemented as single color Sech pulses. 
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FIG. 7: (a): Schematic representation of the controlled phase 
gate implementation of Table ITvl without refocussing pulses. 
Pulses 1-6 implement a controlled phase gate and pulses 13- 
18 implement a 'controlled identity' operation which are ro- 
bust to variations in A*. Taken together pulses 1-6 and 13-18 
are robust to variations in both A c and A t . 
(b): Representation of the controlled phase gate sequence 
with refocussing. Note that since the dipole blockade pro- 
hibits simultaneous excitation of two interacting qubits the 
NOT gates implemented by pulses 7-12 and 19-24 must be 
staggered with respect to one another. 



the bit-flip approach to rephasing to work best when the 
qubit is only mildly perturbed between each flip. Since 
almost all non-refocussing manipulations in the gate se- 
quence address the target qubit it is not surprising that 
it is more vulnerable to the effects of the hyperfine shift. 
If the simulation is run without the refocussing pulses 
one obtains fidelities over the range T =0.1-1. 



V. DECOHERENCE 

For single pulses and short pulse sequences the effects 
of decoherence can usually safely be ignored. The pro- 
posed pulse sequence of Table IIVI is relatively long and 
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FIG. 8: Fidelity of the controlled phase gate implementation 
against the hyperfine shift of the control (8 C ) and target (St) 
qubits. The initial state is chosen to be the equally weighted 
superposition state \ip in ) = i (|00) + |01) + 1 10) + |11>) and 
the pulse parameters are flo = 4 MHz, (3 = 1.28 MHz, p = 
3, T = 1.5 /is, A c = 100 kHz and A< = 80 kHz. Note the 
sequence is virtually insensitive to variations in S c . 



FIG. 9: Fidelity of the controlled phase gate implementation 
including decoherence effects. The initial state is chosen to 
be the equally weighted superposition state \tpi n ) = §(100) + 
01) + 1 10) + 1 1 1 ) ) and the pulse parameters are fio = 4 MHz, 
(3 = 1.28 MHz, p = 3, T = 1.5 ps, A c = 100 kHz and At = 80 
kHz. 



for this reason it is worthwhile to consider decoherence 
effects in simulations. In this section we identify the main 
sources of decoherence in the rare earth quantum com- 
puting setup. The hyperfine ground states have very long 
coherence times, as long as 82 ms has been reported [ll| . 
The main source of decoherence is thus attributable to 
populating the excited state. 

In order to handle the dynamics in the face of relax- 
ations we turn to the density operator formalism and seek 
to solve the master equation 

p = i[p,H] + C ielax (p). (22) 

The function Cr S ^(p) is required to be on the so called 
Lindblad form [T^l 

Aeto(p) = -f Y,{ClCmP + P ClC m } + ]T C m pCj,. 

m m 

(23) 

The operators C m describe the relevant decoherence pro- 
cesses. We shall assume that the dominant source of de- 
coherence is spontaneous decay from the excited state. 
We thus include the following operators for each ion 

C { spm) = v^T|0><e| 

c (s P on) = v /^f (24) 

The two operators describe spontaneous transitions from 
the excited state to the two qubit states with branching 
ratios bo and b\. 

Finally we solve numerically Eq. H22JI for the full con- 
trolled phase gate pulse sequence. A plot of the fidelity 
as a function of the hyperfine shifts of the control and 
target qubits is given in Fig. We observe a drop in the 
fidelity of the gate sequence compared to Fig. |SJ How- 
ever, we believe we have correctly identified and included 



the prominent source of decoherence and so the fidelities 
of Fig. El should closely resemble those obtained in the 
laboratory. In principle the fidelity could be further in- 
creased by increasing the laser intensity thus allowing 
faster pulses while remaining in the adiabatic limit. This 
would be pushing the limits of current technical feasibil- 
ity so we do not explore this avenue further. 

It is interesting to analyze the connection between the 
fidelity and the lifetime T e . In Fig. ^] we plot the gate fi- 
delity as a function of the hyperfine shift of the target ion 
for different lifetimes (recall from Figs. 03 and [5] that the 
fidelity is all but independent of the control qubit hyper- 
fine shift) . We observe a marked improvement in fidelity 
as T e becomes longer, a point which emphasizes the im- 
portance of choosing an ion with the right properties for 
the experimental realisation. 

VI. CONCLUSION 

In conclusion we have investigated the effects of ran- 
dom but constant shifts of the hyperfine qubit levels in 
the REQC system. The evolution of a two-level system 
during a complex hyperbolic secant pulse has been anal- 
ysed analytically to first order using perturbation theory. 
We have used these results to explore the consequences of 
the hyperfine shift during single qubit rotations and have 
found that errors are dynamically suppressed due to the 
AC Stark shift. We have presented a revised implemen- 
tation of the controlled phase gate which is resilient to 
the hyperfine shift and have shown that decent fidelities 
are possible even when considering the effects of decay 
and decoherence. 

The fidelities are not quite as high as could be desired 
however this is chiefly due to relaxation and decoherence 
effects caused by prolonged integrated population of the 
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FIG. 10: (Color online) Fidelity as a function of the target ion 
hyperfine shift for different excited state lifetimes T e . Starting 
from below the solid curve is for T e = 100 /xs, the dashed curve 
for T e = 500 /j,s, the dotted curve for T e = 1 ms and the bold 
curve for T e = 1 s. 



excited state which is inevitable if we are to correct for 
the errors inherent in the system. Recently Wesenberg et 
al. have proposed a single instance approach to REQC 
|13| which would eliminate some of these errors thereby 
allowing higher fidelities. 

The gate implementation we have proposed is con- 
structed from intuitive primitives and there is a clear 
physical idea behind each pulse. Presumably a more 
efficient implementation could be found using optimal 
control theory, possibly using our sequence as a starting 
point for numerical optimisation. 
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